Preprint typeset using I4TgX style emulateapj v. 5/2/1 1 



LATE-TIME LIGHT CURVES OF TYPE II SUPERNOVAE: PHYSICAL PROPERTIES OF SNE AND THEIR 

ENVIRONMENT 

Masaaki Otsuka 1,2 , Margaret Meixner 1 , Nino Panagia 1 , Joanna Fabbri 3 , Michael J. Barlow 3 , Geoffrey C. Clayton 4 , 
Joseph S. Gallagher 5 , Bene. K. Sugerman 6 , Roger Wesson 3 , Jennifer E. Andrews 4 , Barbara Ercolano 7 , and 

Douglas Welch 8 
(Received; Revised; Accepted) 

ABSTRACT 

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 
2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope {HST), Gem- 
ini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up 
to day 575 because it was serendipitously in our SN 2002hh field of view. Of our 43 HST observations in 
total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of 
all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and 
around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as 
well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate 
decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data 
indicate that dust is forming in our targets from day ^300-400, supporting SN dust formation theory. In this 
paper we focus on other physical properties derived from the late time light curves. We estimate 56 Ni masses 
for our targets (0.5-14 x 10" 2 M©) from the bolometric lightcurve of each for days ^150-300 using SN 1987A 
as a standard (7.5xl0~ 2 M ). The flattening or sometimes increasing fluxes in the late time light curves of 
SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar 
hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively 
dense materials (n(H) >400 cm" 3 ) and SNe 2003gd and 2004et have densities more typical of the interstellar 
medium (~1 cm" 3 ). We analyze the sample as a whole in the context of physical properties derived in prior 
work. The 56 Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10" 2 , supporting the 
previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear 
to be correlated with progenitor mass. 

Subject headings: (Stars:) supernovae: individual: SN1999bw, SN2002hh, SN 2003gd, SN2004et, SN2005cs, 
SN2006bc, SN 2008S; (ISM:) dust, extinction 



1. INTRODUCTION 

Since the discovery of SN 1987 A, dust formation has 
been confirmed observationally in the ejecta of several core- 
collapse supernovae (Type II SNe; e.g., Wooden et al. 1993, 
Kozasa et at al. 2009). In these, one or more of the follow- 
ing signatures are observed ^250-600 days after the initial 
explosion: (1) the development of an infrared excess seen in 
the spectral energy distribution (SED) of the SN; (2) the ap- 
pearance of blue-shifted optical/near-infrared lines, which is 
interpreted as the attenuation of the redshifted wing (i.e., the 
receding gas) by dust formed homogeneously in the ejecta; 
and (3) the steepening decline of the optical lightcurve con- 
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sistent with an increase in circumstellar extinction. The SN 
luminosity after day ~ 150 is supported by 7 rays from 56 Co 
decaying into 56 Ni. If the dust is formed in SNe, the opti- 
cal luminosity from the SN remnant has become much fainter 
than would be produced by the radioactive decay power. 

In order to properly assess dust formation by Type II SNe 
ejecta it is essential to continue monitoring Type II SNe be- 
fore and after dust formation started (i.e. day ^300) because 
other physical properties such as the 56 Ni mass and light echos 
from circumstellar and interstellar dust can be measured and 
disentangled from dust production effects. The early epoch 
evolution of Type II SNe before day ^300 has been relatively 
well studied in optical and near-IR bands. However at the 
present time, the evolution beyond day 300 has been poorly 
measured. 

To address these issues, we have conducted a long term 
photometric monitoring of SNe using high angular resolution 
facilities such as HST/WFPC2/NICMOS, WIYN/WHIRC, 
and Gemini/NIRI through our SNe project SEEDS (Survey 
for Evolution of Emission from Dust in SNe, PI: M. Barlow). 
In this paper we report on the late-epoch light-curves of 6 
Type II SNe. The remainder of this paper is organized as fol- 
lows. In section 2, we present the observations, data process- 
ing and resulting images, and light curves. In section 3, we 
analyze the resulting light curves to derive 56 Ni masses for 
all targets except for SN 1999bw and, for cases with detected 
light echoes, to estimate the density of the environment. In 
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section 4, we discuss the results for each SN. In section 5, we 
do a comparative analysis of these SNe in the context of phys- 
ical properties reported in the literature, such as the amount of 
dust produced and progenitor mass. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. Sample selection 

We chose 6 Type II SNe in nearby galaxies, matching the 
following three criteria. (1) The distance to a SN is within 
^15 Mpc and thus suitable for long-term follow-up observa- 
tions with HST, Spitzer, and >4-m class ground based tele- 
scopes. (2) At the start of the observations, the SN is 1 year 
or older after the initial explosion. (3) The expected V-band 
magnitude of targets is <28 based on the evolution of SN 
1987 A in order that the SNe is detectable within a reason- 
able integration time. These 6 Type II SNe include 4 Type 
II-P (SNe 2002hh, 2003gd, 2004et, and 2005cs), 1 Type Il-n 
(SN 1999bw), and 1 Type II-L or II-P (SN 2006bc). In addi- 
tion, during our campaign, a supernova impostor, SN 2008S, 
erupted in the field of view of our SN 2002hh observations, 
allowing us to monitor this new type of variable as well. 

2.2. HST observations 

The high-angular resolution images of our target SNe were 
taken by the Hubble Space Telescope with the Wide Field 
Planetary Camera 2 (WFPC2) and the Near Infrared Camera 
and Multi-Object Spectrometer 2 (NICMOS2) through the cy- 
cle 16 proposal by M.Meixner (proposal ID: GO 1 1229). The 
observations were coordinated so that the WFPC2 and NIC- 
MOS images were taken within days to a week of each other. 
In the WFPC2 observations, we took the images using the 
F450W, F606W, F622W, and F8 14W filters. In the NICMOS2 
observations, we took the images using the F110W, F160W, 
and F205W filters. The observation logs and SN positions are 
summarized in Table Q] The boldface coordinates in this ta- 
ble are derived from our F606W or F622W HST observations 
and the others are from the SIMBAD database. 

2.3. HST data reduction 

The data reduction and calibration were performed by stan- 
dard techniques. Using the STSDAS package version 3.8 
specialized for HST, we reduced and calibrated the data, in- 
cluding cosmic rays and artificial features removal and de- 
tector linearity correction. The final high resolution images 
were created using the STSDAS/DRIZZLE package. For both 
the WFPC2 and NICMOS2 observations, we conducted sub 
pixel scale dithering to improve S/N and remove cosmic rays 
and obtain further high-resolution images via the drizzle tech- 
nique. We set 4 or 5 points dithering with 115-400sec per 
science frame for each band. After DRIZZLE processing, 
we performed an additional distortion correction and align- 
ment using background stars for the NICMOS2 data. The 
point-spread-function (PSF) of these reference background 
stars were fit by Gaussian profiles to obtain accurate posi- 
tions, and then instrumental distortions were corrected with 
the XYXYMATCH, GEOMAP, and GEOTRAN IRAF rou- 
tines, and we aligned the positions of stars in the NICMOS 
detector to those in the WFPC2. The resultant plate scale 
reaches 0.025" pixel" 1 in the WFPC2 and archived ACS/HRC 
data (see below) and 0.04" pixel" 1 in the NICMOS2 and 
archival NICMOS3 data (see below), respectively. The size of 
the PSF at FWHM typically corresponds to ^3 pixel. Using 
the PSF fitting package IRAF/DAOPHOT, photometry was 



performed on the images, except for SN 2005cs in F110W 
and F160W and SN 2002hh in all bands. For SN 2005cs 
in Fl 10W and F160W, we used differential images with pre- 
explosion images using archival NICMOS3 images of the host 
galaxy taken in 1998 June 28th (PLN.Scoville, Proposal-ID: 
GO 7237) and our NICMOS2 data to detect the SN alone and 
performed PSF fitted photometry. SN 2002hh has an extended 
shell structure as shown in Fig.l, we therefore performed 
aperture photometry based on background star subtracted im- 
ages. We used 0.6" radius regions to define the entire the SN. 
We subtracted background from an annulus centered on the 
SN with inner and outer radii of 0.6" and 0.8". The aperture 
corrections were calculated with an empirical PSF function. 
Note that the magnitudes of SN2002hh listed in Tables 1, 2, 
and 3 include the SN + extend light echo. 

For each filter, the measured count rate (CR, in units of 
DN s" 1 ) of the SN was converted to flux by multiplication 
with the PHOTFLAM (erg s" 1 cm" 2 A" 1 DN" 1 or erg s" 1 
cm" 2 A" 1 electron" 1 ) conversion factor given in the fits header, 
where PHOTFLAM is the bandpass-averaged flux density for 
a source that would produce a count rate (CR in units of DN 
s" 1 or electron s" 1 ). The HST system magnitude (STMAG) 
was calculated from the following equation, 

STMAG = -2.5 log 10 (C# x PHOTFLAM) + ZP, ( 1 ) 

where ZP is the zero point magnitude for 1 DN s" 1 or 1 elec- 
tron s" 1 . For the WFPC2 and archival ACS data (see be- 
low), the conversion from F450W, F606W and F622W, and 
F814W STMAG to the Vega magnitude system B, V, R, and 
/, respectively, used the STSDAS/SYNPHOT package. The 
broad band NICMOS filters, F110W, F160W and F205W, 
are roughly equivalent to the 7, H and K filters respectively, 
color transformations are not well constrained for late-time 
SN spectra at near-IR wavelengths. Based on our experi- 
ence, for SNe in days 300-1000, the spectral energy distri- 
bution (SED) in the range from ^0.44 to ^2.2 /im can be 
represented by a single blackbody function with the effec- 
tive temperature of ^7 000- 1 1 000 K. The conversion from 
the STMAG to the Vega magnitude system was done by the 
STSDAS/SYNPHOT package, assuming a 9 000 K blackbody 
function as the incident SED. 

2.4. HST Archive data 

Except for SN 2004et, archival WFPC2, Advanced Cam- 
era for Surveys (ACS), and NICMOS3 images are available 
for our SNe sample as noted in Table 2. We created drizzled 
images and we measured the magnitudes by the same way as 
described above. The measured magnitudes are listed in Table 
12 In this table we also list the magnitudes measured at earlier 
epochs as found in the literature as noted in the table. 

2.5. NIRI and WHIRC observations 

To complement HST observations, we performed JHK- 
band high-resolution imaging observations for SNe 2002hh, 
2004et, and 2008S using the Gemini north 8.2-m telescope 
with the Near Infrared Imager and spectrometer (NIRI) and 
the WIYN 3.5-m telescope with the WIYN High Resolu- 
tion Infrared Camera (WHIRC; Meixner et al. 2010). The 
NIRI and Aug 2009 WHIRC observations were carried out 
through the NOAO open use programs of Pis. G.C.Clayton 
and M. Otsuka, respectively. The other WHIRC observations 
were carried out in STScI guaranteed observation time (PL: 
M.Meixner). The plate scale of both instruments is ^0.1" 
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TABLE 1 

HST Observation log and measured magnitudes 



SN Name 


Host Galaxy 


D(Mpc) 


E(B-V) 


q(2000.0) 


(5(2000.0) 


Obs.Date 


Epoch (day) 


Exp. time 


Detector 


Filter 


Detect 


Magnitudes 


SN1999bw 


NGC3198 


13.7 


0.013 


10:19:46.81 


+45:31:35.0 


2008/04/03 


3271.24 


128secx5 


NIC2 


F110W 


N 




(ll-n) 








10:19:46.68 


+45:31:37.32 


2008/04/03 


3271.25 


128secx4 


NIC2 


F160W 


N 
















2008/04/03 


3271.26 


115secx5 


NIC2 


F205W 


N 
















2008/04/03 


3271.11 


400sec x 4 


PCI 


F606W 


Y 


25.77 ± 0.16 














2008/04/03 


3271.20 


400sec x 4 


PCI 


F814W 


Y 


25.55 ± 0.24 


SN2002hh 


NGC 6946 


5.9 


1.97 


20:34:44.29 


+60:07:19.0 


2007/07/09 


1716.88 


400sec x 4 


PCI 


F606W 


Y 


20.62 ± 0.11 


(II-P) 








20:34:44.19 


+60:07:20.11 


2007/07/09 


1716.91 


400sec x 4 


PCI 


F814W 


Y 


19.67 ± 0.11 














2007/07/09 


1716.96 


400sec x 4 


PCI 


F450W 


Y 


23.25 ± 0.18 














2007/07/10 


1717.02 


128secx5 


NIC2 


F110W 


Y 


18.84 ± 0.11 














2007/07/10 


1717.03 


128secx4 


NIC2 


F160W 


Y 


18.53 ± 0.11 














2007/07/10 


1717.04 


115secx5 


NIC2 


F205W 


Y 


17.44 ± 0.11 


SN20O3gd 


NGC 628 


7.2 


0.18 


01:36:42.65 


+15:44:19.9 


2007/08/11 


1519.79 


400sec x 4 


PCI 


F622W 


Y 


24.91 ± 0.18 


(n-P) 








01:36:42.71 


+15:44:20.25 


2007/08/1 1 


1519.86 


400sec x 4 


PCI 


F814W 


Y 


25.60 ± 0.19 














2007/08/14 


1523.25 


128secx5 


NIC2 


F110W 


Y 


<24.92 














2007/08/14 


1523.26 


128secx4 


NIC2 


F160W 


Y 


<25.47 














2007/08/14 


1523.27 


115secx5 


NIC2 


F205W 


N 




SN2004et 


NGC 6946 


5.9 


0.41 


20:35:25.33 


+60:07:17.7 


2007/07/08 


1071.74 


400sec x 4 


PCI 


F606W 


Y 


23.01 ± 0.12 


(n-P) 








20:35:25.36 


+60:07:17.80 


2007/07/08 


1071.82 


400sec x 4 


PCI 


F814W 


Y 


23.04 ± 0.14 














2007/07/08 


1071.91 


128secx5 


NIC2 


F110W 


Y 


21.76 ± 0.14 














2007/07/08 


1071.94 


128secx4 


NIC2 


F160W 


Y 


22.56 ± 0.22 














2007/07/08 


1071.95 


115secx5 


NIC2 


F205W 


Y 


<22.07 














2008/01/19 


1266.71 


400sec x 4 


PCI 


F606W 


Y 


23.35 ± 0.11 














2008/01/19 


1266.78 


400sec x 4 


PCI 


F814W 


Y 


23.35 ± 0.13 














2008/01/20 


1267.58 


128secx5 


NIC2 


F110W 


Y 


21.99 ± 0.15 














2008/01/20 


1267.59 


128secx4 


NIC2 


F160W 


Y 


22.62 ± 0.24 














2008/01/20 


1267.59 


115secx5 


NIC2 


F205W 


Y 


21.07 ± 0.19 


SN2005cs 


NGC 5194 


7.1 


0.14 


13:29:53.37 


+47:10:28.2 


2007/10/28 


851.96 


400sec x 4 


PCI 


F606W 


Y 


25.16 ± 0.19 


(n-P) 








13:29:52.80 


+47:10:35.64 


2007/10/28 


852.02 


400sec x 4 


PCI 


F814W 


Y 


25.07 ± 0.24 














2007/10/31 


854.15 


128secx5 


NIC2 


F110W 


Y 


<23.30 














2007/10/31 


854.16 


128secx4 


NIC2 


F160W 


N 
















2007/10/31 


854.17 


115secx5 


NIC2 


F205W 


N 
















2008/03/31 


1006.15 


400sec x 4 


PCI 


F606W 


Y 


<25.94 














2008/03/31 


1006.21 


400sec x 4 


PCI 


F814W 


N 




SN2006bc 


NGC 2397 


16 


0.52 


07:21:16.50 


-68:59:57.3 


2007/09/07 


531.98 


400sec x 4 


PCI 


F606W 


Y 


23.16 ± 0.19 


(n-P or L) 








07:21:16.45 


-68:59:56.75 


2007/09/07 


532.07 


400sec x 4 


PCI 


F814W 


Y 


21.72 ± 0.19 














2007/09/17 


541.45 


128secx5 


NIC2 


F110W 


Y 


21.60 ± 0.20 














2007/09/17 


541.46 


128secx4 


NIC2 


F160W 


Y 


20.80 ± 0.20 














2007/09/17 


541.47 


115secx5 


NIC2 


F205W 


Y 


19.90 ± 0.10 














2008/02/17 


694.25 


400sec x 4 


PCI 


F606W 


Y 


23.85 ± 0.26 














2008/02/17 


694.28 


400sec x 4 


PCI 


F814W 


Y 


22.51 ± 0.31 














2008/02/17 


694.35 


128secx5 


NIC2 


F110W 


Y 


22.40 ± 0.30 














2008/02/17 


694.36 


128secx4 


NIC2 


F160W 


Y 


21.70 ± 0.40 














2008/02/17 


694.37 


115secx5 


NIC2 


F205W 


Y 


20.70 ± 0.20 



REFERENCES. — (For D and E(B-V)) SN 1999bw: Freedman et al. (2001), Smith et al. (2010); SN 20()2hh: Barlow et al. (2005), Karaehentsev et al. (2000), Pozzo et al. (2006); 
SN 2003gd: Van Dyk et al. (2003), Hendry et al. (2005); SN 2004et: Fabbri et al. (201 la), Karaehentsev et al. (2000); SN 2005cs: Takats, K., & Vinko (2006), Tsvetkov et al. (2006); 
SN 2006bc: Gallagher et al. (2010); Mihos & Bothun (1997) 

NOTE. — The measured magnitudes at F450W, F606W, F622W, F814W, Fl 10W, F160W, and F205W in column (10) are converted to B, V, R, I, J, H, K magnitudes. We converted 
the STMAG into the Johnson-Cousins system and Vega system using the STSDAS/SYNPHOT In the fifth and sixth columns: the boldface characters are the SN positions derived from 
our HST observations and the other characters are from SIMBAD database. In the tenth column, PCI and NIC2 mean the PCI chip in the WFPC2 and the NICMOS2 chip in the 
NICMOS, respectively (see text for detail). 



pixel -1 and the seeing was 0.3"-0.8" at JHK-bands. The av- 
erage PSF size of ^0.5" enabled us to measure magnitudes 
of both SNe with relatively little contamination from nearby 
stars. 

The data reductions using the IRAF/MSCRED package 
were carried out in a standard manner (dark subtraction, flat- 
fielding, and local sky subtraction). For the WHIRC data, we 
performed a detector linearity correction using the WHIRC 
task WPREP and distortion correction using the files down- 
loaded from the WIYN/WHIRC web page. The magnitude 
measurements were performed using the IRAF/DAOPHOT. 
The conversion of the instrumental magnitudes into the stan- 
dard system was performed relative to Two Micron All-Sky 
Survey (2MASS) photometry of 4-10 stars in the field. 

The measured magnitudes are listed in Table [3] The er- 
rors include the estimated errors of 2MASS stars' magnitudes. 



For SN 2002hh, we converted the Vega magnitudes derived 
from HST data to 2MASS magnitudes. Comparison between 
HST and 2MASS magnitudes indicates that a star nearby SN 
2002hh contributes ^0.2-0.3 mag. 

2.6. Results 

The resultant images of target SNe are presented in Fig[T] 
We are successful in detecting the SNe and measuring either 
F606W(V) or F622W(fl) band magnitudes of all the targets 
in HST observations. The light-curves of each SN are pre- 
sented in Figsf2] (VRI) and |3] (JHK-bands). Theses light- 
curves are based on our data (filled circles) and data from 
the literature for each object (open circles). The data after 
>600-1000 days are mainly from the former and the data be- 
fore this phase are mainly from the latter. The data sources for 
each object are summarized in Table|4] The values in Figs. [2] 
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FIG. 1.— The HST images of SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and2006bc. The SN is in the center of each panel. The size of each panel is l"x 
l"and each tick mark is 0.2". The intensity scale (erg s cirT 2 A -1 pixel -1 ) is common in each SN. In all cases north is up and east is to the left. Each box has 
a label for HST filter and day after explosion. 



and |3] are de-reddened absolute magnitudes adopting the dis- 
tance to the host galaxy, D, and E(B-V) listed in the third 
and forth columns of Table 1. All E(B-V) in Table 1 are 
the total extinction by the Milky Way Galaxy and host galaxy. 
We used the reddening function assuming the Ry=3. 1 case by 
Cardelli et al. (1989) for all SNe. As references, we plot the 
lightcurves of the Type II-P SNe 1987A (Hamuy & Suntzeff 
1990; Walker & Suntzeff 1991) and 1999em (Elmhamdi et al. 
2003), which are well studied Type II SNe. 

We measured the rate of fading of SNe jband during 4 
phases, namely, days ^30-100 (plateau phase), ~ 150-300 (ra- 
dioactive decay phase), ~300-800 (dust forming phase), and 
>800. The results are summarized in Table[5] For each phase, 
we fit a line to the light-curve, the slope of which is the fad- 
ing rate for that phase measured in units of magnitude per 100 
days. The values indicated by bold face are derived from the 
data presented in this work. We note that light echoes, which 
may appear as flattening or increasing fluxes in the lightcurve, 



could arise during any of these phases depending on the loca- 
tion of the scattering material with respect to the SNe. The 
lightcurves of Type II SNe during the radioactive decay phase 
can be explained by 7-rays power by the radioactive decay of 
56 Co to 56 Fe, at a rate corresponding to the e-folding life time 
of the 56 Co decay (t^ = 111.3 days). The rate expected for 
radioactive decay of 56 Co into 56 Fe is 0.98 mag per 100 days, 
while the averaged 7^ among 8 SNe is 0.96±0.12 (Patat et 
al. 1994). For SN 2004et, for example, Sahu et al. (2006) 
found that the BVRI lightcurves in 180-310 days were linear 
with 7^ = 1.04, 7r = 1.01, 7/ = 1.07. Maguire et al. (2010) 
also obtained similar results (1.02, 0.92, 1.09, respectively). 
Our estimated jvri except for SNe 2002hh and 2005cs are in 
good agreements with Patat et al. (1994). Two SNe might 
have additional power sources in this phase (See sections 4.2 
for SN 2002hh and 4.5 for SN 2005cs). We find that y HK are 
~1-1 .4 for all objects. Beyond ^300 days, the fading rates 
are steeper than the 56 Co decay, because the opacity to 7-rays 
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TABLE 2 

The magnitudes by previous HST observations. 



SN Name 


Obs.Date 


Epoch 


Detector 


Filter 


Magnitudes 


Ref. 


SN 1999bw 


2001/01/09 


629.88 


PCI 


F555W(V) 


24.08 ± 0.06 


(1) 


SN 2002hh 


2005/09/17 


1056.26 


HRC 


F606W(V) 


20.88 ±0.11 


(2) 




2005/09/17 


1056.40 


HRC 


F814W(/) 


19.49 ±0.11 


(2) 




2006/04/23 


1273.94 


HRC 


F606W(V) 


20.91 ±0.11 


(2) 


SN 2003gd 


2004/12/08 


544.11 


HRC 


F435W05) 


23.98 ±0.15 


(3) 




2004/12/08 


544.11 


HRC 


F622W(fi) 


23.13 ±0.16 


(3) 




2007/06/21 


1469.45 


PCI 


F450W(B) 


25.01 ±0.17 


(4) 




2007/06/21 


1469.45 


PCI 


F675WCR) 


24.10 ±0.16 


(4) 


SN 2005cs 


2006/12/10 


529.02 


HRC 


F555W(V) 


22.07 ±0.15 


(5) 




2006/12/10 


529.02 


HRC 


F814W(/) 


21.36 ±0.15 


(5) 


SN 2006bc 


2006/10/14 


203.28 


WFC1 


F555W(V) 


22.10 ±0.10 


(6) 




2006/10/14 


203.28 


WFC1 


F814W(/) 


19.80 ±0.10 


(6) 



NOTE. — HRC means High Resolution Camera mode in the ACS. WFC1 means the 
wide field camera 1 chip in the WFPC2. 

REFERENCES. — (1) Li et al. (2002); (2) This work (data source; GO 10607, PI: 
Sugerman); (3) Sugerman (2005); (4) Maund & Smartt (2009); (5) This work (data 
source; GO 10877, PI: Li); (6) Gallagher et al. (2011; GO 14980, PI: Smartt) 



TABLE 3 
NIRI AND WHIRC OBSERVATIONS. 



SN 2002hh 



Obs.Date 


Epoch 


Instr. 


J 


H 


K 


2005/06/04 


951.5 


NIRI 


18.50 ±0.16 


17.76 ±0.16 


16.93 ± 0.20 


2005/08/04 


1012.5 


NIRI 


18.57 ± 0.15 


17.80 ±0.16 


16.93 ± 0.20 


2005/10/16 


1085.3 


NIRI 


>18.48 


> 17.88 


16.78 ± 0.20 


2006/07/07 


1349.5 


NIRI 


18.51 ± 0.15 


17.75 ± 0.16 


16.95 ± 0.20 


2007/07/10 


1717.0 


HST 


18.61 ± 0.19 


18.03 ±0.23 


17.02 ± 0.21 


2008/11/12 


2208.1 


WHIRC 


18.59 ±0.15 


17.82 ±0.16 


17.00 ± 0.20 


2009/05/10 


2387.4 


WHIRC 




17.54 ±0.20 




2009/08/30 


2499.3 


WHIRC 




18.07 ±0.20 










SN 2004et 






2009/08/30 


1801 


WHIRC 




<22.6 










SN 2008S 






2008/11/12 


283.3 


WHIRC 


18.37 ± 0.15 


16.73 ±0.15 


15.64 ± 0.20 


2009/05/10 


463.5 


WHIRC 




17.83 ±0.16 


16.25 ± 0.73 


2009/08/30 


575.5 


WHIRC 




18.74 ±0.16 





NOTE. — The errors 
the HST magnitudes wer 
magnitude. 



of 2MASS stars' magnitudes are included. In SN 2002hh, 
e re-calibrated using the NIRI data and nearby 2MASS stars' 



TABLE 4 

The data source of Type ii SN and SN 2008S lightcurves. 



SN Name 



References 



SN 1999bw Li et al. (2002); Smith et al. (2010); This work 

SN 2002hh Pozzo et al. (2006); Welch et al. (2007); This work 

SN 2003gd Hendry et al. (2005); Maund & Smartt (2009); This work 

SN 2004et Fabbri et al. (201 la); Kotak et al. (2009); This work 

SN 2005cs Pastorello et al. (2009); Tsvetkov et al. (2006); This work 

SN 2006bc Brown et al. (2009); Gallagher et al. (2010); This work 

SN 1987A Hamuy & Suntzeff (1990); Walker & Suntzeff (1991) 

SN 1999em Elmhamdi et al.(2003) 

SN 2008S Botticella et al. (2009);Prieto et al. (2010); This work 

is decreased and dust starts forming. This supports SN dust 
formation theory. We find flattening or increasing magnitude 
features in SNe 1999bw, 2002hh, 2003gd, 2004et, and 2006bc 
at the same epoch. We will discuss whether these SNe have 
light echos based on the surrounding ISM density later. 

3. ANALYSIS OF LATE TIME LIGHT CURVES 

In our prior work on the late time light curves of these SNe, 
we have concentrated on the dust formation process and es- 
timated dust mass as noted in section 4. In this paper, we 
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FIG. 2. — VRI lightcurves on a de-reddened absolute magnitude scale. The 
open circles are from published papers and the filled circles indicate our con- 
tributions. The broken lines are from published SN 1987A data. The data 
sources are listed in Tablef4] 

analyze these light curves for two other quantities: the ejected 
56 Ni mass (Tabled and the density of the environment when 
a light echo is present (Table [8j. Table [8] summarizes the de- 
rived physical parameters for these SNe including 56 Ni mass, 
the environment density, dust produced and progenitor mass. 
Below we describe how we calculate the 56 Ni mass and the 
environment density. 

3.1. Mass of 56 Ni 

Here, we focus on the 56 Ni mass synthesized in the SN ex- 
plosion, because we can estimate this mass using bolometric 
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TABLE 5 

The magnitude decline rate (■yband) of Type ii SNe. 



~]band 


Epoch 


SN1999bw 


SN2002hh 


SN2003gd 


SN2004et 


SN2005cs 


SN2006bc 


SN1987A 


SN1999em 


TV 


30-100 








0.764 


0.067 






0.400 




150-300 




0.972 


0.780 


1.006 


0.277 


0.709 


0.980 


0.932 




300-800 








1.268 


1.069 


0.333 


1.469 


0.985 




>800 


0.064 


0.008 




0.173 


0.506 




0.562 




JR 


30-100 




0.147 




0.417 


-0.591 






-0.030 




150-300 




0.971 


1.211 


1.014 


0.684 




0.855 


0.945 




300-800 






1.453 


1.509 


1.135 




1.640 


1.126 




>800 






0.052 








0.528 




7i 


30-100 






0.798 


0.105 


-0.856 






-0.269 




150-300 




0.968 


1.295 


0.918 


0.771 


0.669 


0.946 






300-800 




0.135 




1.454 


1.002 


0.468 


1.661 


1.126 




>800 




0.049 


0.052 


0.126 






0.565 




■yj 


30-100 








-0.467 


-1.164 






-0.371 




150-300 




1.256 




1.497 


1.071 




1.192 


1.158 




300-800 








1.016 


1.017 


0.493 


1.528 






>800 




0.007 




0.117 






0.809 




7H 


30-100 
150-300 
300-800 
>800 




1.316 
0.014 




-0.269 
1.438 
1.098 
0.003 


-0.911 


0.555 


1.305 
1.386 
0.857 


-0.428 
1.444 


7K 


30-100 
150-300 
300-800 
>800 




1.286 
0.011 




-0.245 

1.177 
-0.510 


-0.784 
1.790 


0.493 


1.578 
1.594 
0.624 


-0.425 
1.482 



Note 

observations and the archival HST data 



The values indicated by bold face are derived from the data taken by our 



lightcurves during radioactive decay. 

The method for 56 Ni mass estimation is to (1) determine the 
bolometric luminosity light curve for the SNe and (2) com- 
pare the slope during radioactive decay between targets and a 
standard SN. Here, we use SN 1987 A as a standard. Arnett & 
Fu (1989) estimated 0.073 ± 0.015 M of freshly synthesized 
56 Ni in SN 1987 A by comparing the observed bolometric lu- 
minosity curve and theoretical model. 

By using simple blackbody fitting techniques, we estimated 
the luminosity Lhot and emitting radius /?hot of the hot compo- 
nent r e ff (>5000 K here). We estimated these three quan- 
tities based on BVRIJHK. Actually, there are archival UV 
data only in very early phase. For SN 2004et, Fabbri et al. 
(201 1) performed multi components blackbody fittings to the 
BVRIJHK and Spitzer IRAC 4 bands and MIPS 24/im data. 
They showed that the percentages of hot component to the 
total luminosity exceeded 76 % from days 64 to 464. We as- 
sume that the condition is the same as in the other SNe, i.e., 
the hot component represented by optical- to near-IR SED is 
the main component of the bolometric luminosity. Since we 
are interested in the late epoch spectral energy distributions 
(SEDs), we selected days 150, 300, and 600. The flux densi- 
ties at each band are estimated using the fading rate in Table 
[5] The results of SED fittings are summarized in Table|6] The 
uncertainties of Lhot and T e ff are ^20% of the estimated values 
and ± ~500 K. 

Next, we estimated the ejected 56 Ni mass based on the as- 
sumption that the energy source of our sample SNe beyond 
150 days is mainly from the power of gamma rays emitted 
through 56 Co to 56 Fe radioactive decay. As the ejecta expand, 
their opacity to 7-rays is expected to decrease, which results 



in a modified light curve given by Woosley et al. (1989), 



Ltotai = Const x exp 



?56 



1 — exp 



-^56, 



^ (7) 



(2) 

where tse is the e-folding life time of 56 Co and the term in the 
brackets is the effective opacity; Kg6, 7 = 0.033 cm 2 g" 1 is the 
average opacity to 56 Co-decay 7-rays, and (j>o = 7 x 10 4 g cm -2 
is the column depth at the fiducial time 1 = 1 1.6 days. Fitting 
equation (|2]i to the estimated bolometric luminosity curve, we 
can determine the constant for each SN. Assuming that ra- 
dioactive materials were powering the late time photometric 
evolution, we estimated the 56 Ni mass by comparing the lu- 
minosity curves of each SN and SN 1987 A by the following 
equation 

M» m {SN) = Mx, m ($lA) x ConstsN M (3) 
Const &1A 

If we can estimate luminosities at one epoch during the ra- 
dioactive decay phase, we can estimate the 56 Ni mass in prin- 
ciple. The estimated 56 Ni mass is listed in the last column of 
Table|6] Since several bands are lacking for SN 2006bc during 
the radioactive decay phase, the estimated 56 Ni mass is based 
only upon the HST data. Since the quantities for SN 2002hh 
are derived from the bolometric luminosity and flux densities 
of the SN and the light echo, the estimated values might have 
>30 % uncertainties. As later discussed, the 56 Ni mass of SN 
2002hh might go down by ~80 % (i.e., ~2.8 M Q ) when we 
consider the contribution from the light echo to the total flux. 

3.2. Light echoes and environment densities 
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FIG. 3. — The lightcurves at JH AT-bands in de-reddened absolute magnitude 
scale. The symbols and lines are the same meanings in Fig. 1. The data 
sources are listed in Tableg] 

Light echoes are created when the flash from the SNe scat- 
ters off dust in their environments, which can be interstellar 
or circumstellar. The timing of the light echo depends on the 
location of the dust with respect to the SN. However, light 
echoes at late times in the SN light curve are more easy to de- 
tect because the increased light due to the SN is more apparent 
when the SN light has decayed significantly. 

The intensity of the light echo depends on the density of 
the surrounding medium. Using the scattering optical depth 
of the light echo given by Romaniello et al. (2005), 



r ec / !0 , JCfl (A) = 7.6 x 10" 



5 5.8 x 10 21 fl(A) St SN (X) 
k 3~X 100 days 



TABLE 6 

The luminosity, temperature, emitting radius, and 56 Ni mass 



SN Name 


Epoch 
(davl 


Lh t 
(X 10 6 Ln) 


(K) 


Rh t 
(xlO 13 cm) 


M( 56 Ni) 
(XlO -2 Mnl 


SN2002hh 


150 


238 


5490 


119 


14 




300 


39.9 


6170 


38.6 




SN2003gd 


300 


2.4 


5930 


10 


0.8 


SN2004et 


150 


92.9 


7040 


45.1 


5.9 




300 


12.6 


7200 


15.9 






600 


0.3 


6360 


3.26 




SN2005cs 


300 


1.50 


7250 


5.41 


0.5 


SN2006bc 


600 


0.3 


6800 


2.84 


2.7 


SN1987A 


150 


86.5 


6030 


59.4 


7.3 




300 


25.8 


7130 


23.2 






600 


0.4 


6500 


3.42 




SN1999em 


150 


22.0 


5560 


35.2 


1.9 




300 


6.1 


5220 


21.1 





Note. 



- The uncertainties of Lw and T e ffare ' 
w(A) nTk) n(b) 



'20% and ~50() K. 



(4) 



0.6 lcm- 3 M(O) ' 

where k = NhIE{B - V) is the ratio of total neutral hydrogen 
column density to color excess at solar metallicity (£=5.8 x 
10 21 cm" 2 mag -1 ; Bohlin et al. 1978), (5fsN is the duration of 
the burst of the SN (—100 days in the case of Type II-P), w(A) 
is the grain albedo (—0.6 at optical wavelength; e.g., Mathis 
et al. 1977). Using equation (2), T ec h ,sca(V) is estimated to 
be 1.205 xlO" 4 n H adopting uj(V)=0.6 and <5f SN (A)=100 days. 
Sparks (1994) argued that the difference between peak (Msn) 
and light echo magnitudes (M ec ho) is 



M, 



echo : 



• M SN + 0.5-2.5 log(T ech0iSCa ). 



(5) 



If a light echo is associated with a SN, M echo (V) - M SN (V) ~ 
10.3 - 2.51og n(H). Using the light curves, we can measure the 
magnitude difference in SNe with light echoes and estimate 
the density n(H). 

4. RESULTS FOR INDIVIDUAL SNE 
4.1. SN 1999bw 

SN 1999bw is a Type Il-n SN in NGC 3198 (<f=13.7 Mpc; 
Freedman et al. 2001) and was discovered on JD. 245 1288.7. 
The position measured from our HST observation is in agree- 
ment with Sugerman et al. (2004) within 0.5". We de- 
tected the light of the SN in F606W and F814W. The opti- 
cal spectra at early epochs are dominated by emission lines 
having full width at half maximum (FWHM) of -1000-2000 
km s" 1 (Garnavich et al. 1999; Filippenko et al. 1999). 
With Spitzer/IRAC flux densities and the blackbody fitting 
method, Sugerman et al. (2004) estimated the emitting ra- 
dius of 1.6 x 10 16 cm and the expansion velocity of 1000 km 
s" 1 in the fifth year after the initial explosion. Dust production 
has not been measured for this SN so far, most likely because 
insufficient data exist for such an estimate. Fabbri (2011b) 
measured flux densities in 4 Spitzer/IRAC bands (days 3316), 
IRS-PU 16 Aim (3331), and MIPS 24 /im (3320), as listed in 
Table [7] Based on her results, we estimated a dust mass of 
9.2 x 10" 5 M Q at a temperature of -300 K at day ~3320 us- 
ing modified blackbody fitting, assuming amorphous carbon 
grains only with 0.1 /im radius. In the fitting, we adopted the 
optical constants of Zubko et al. (1996). The best fit SED is 
presented in Fig|4j 

Smith et al. (2010) reported BVRI magnitudes at the very 
early plateau phase. They suggested that this object might 
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Wavelength (/mi) 

FIG. 4.— The mid-IR SED of SN 1999bw at days 3316 (3.6/4.5/5.8/8.0 
(im) and 3331 (16 fim). The filled circles are observations. The thick line is 
the result of modified blackbody fitting. 



TABLE 7 
SN 1999BW Spitzer observations 



Obs.Date 


Epoch 
(day) 


3.6 /^m 


4.5 (im 


Flux density (/^Jy) 
5.8 fim 8.0 /^m 


16 /^m 


24 p,m 


2008/05/13 


3316 


<10 


4.0±0.5 


21±3 46±8 






2008/05/17 


3320 










<650 


2008/05/28 


3331 








93±16 





NOTE. — The upper limits for the MIPS 24 /^m data were estimated by measuring 
the flux density of the faint point source seen ~6" to the east of the SN position. This 
was found to have an average brightness of 657 /^Jy with a standard deviation of 2.5 % 
over the four epochs, and this approximate value was used as a robust upper limit for the 
SN brightness at each epoch. This table is taken from Fabbri (201 lb). 



be an impostor from statistical study of expansion velocities 
of Type II SNe and luminous blue variables (LBVs) based 
on optical spectra. They estimated a peak absolute /?-band 
magnitude of -12.65 using optical spectra at day 4, which 
lies between those of P-Cygni and -q Carinae. There is no in- 
formation on the progenitor star. There is no information on 
the brightness of the SN during days 150-626, so we estimate 
the V-band magnitude decline rate for >800 days. Since the 
magnitude difference between peak and day >1000 is —6.6, 
the hydrogen density n(H) assuming solar O abundance (8.69; 
Lodders et al. 2003) is estimated to be —30 cm" 3 , which is 
larger than typical interstellar gas densities (1-10 cm" 3 ). We 
might therefore be looking mostly at pre-existing circumstel- 
lar material rather than the ISM. 

4.2. SN2002hh 

SN 2002hh is a Type II-P SN in a spiral arm of the starburst 
galaxy NGC 6946 (d=5.9 Mpc; Karachentsev et al. 2000) and 
was discovered on JD. 2452574. 5. We detected the light from 
the SN in HST, NIRI, and WHIRC observations. The SN 
birth rate in this galaxy is very high, —0.1 SN yr" 1 during the 
past 100 yr. Within 2' radius from SN 2002hh, SNe 1969P 
and 2008S were observed (Fig(5]). The high SN formation 
rate could be a natural consequence of the high star forming 
activity in NGC 6946. 

SN 2002hh is still bright and close to peak brightness even 
after —7 years (Fig. 3). While at the early plateau phase 
around day 150 the fading rate roughly followed the 56 Co de- 
cay power (0.97 at V and / and — 1 .3 at JHK), the fading rates 



are <0.1 from day 300 or later. Pozzo et al. (2006) argued 
that the near-IR excess in this SN is due to an IR light echo 
from a pre-existing, dusty circumstellar medium. In FigO we 
show the temporal evolution of the SN at F606W (V-band). 
The 2005 and 2006 images (upper left two panels) were taken 
by the ACS/HRC. In the middle three panels, we show the 
subtraction images between two different epochs. Two arc 
structures are seen and appear to be expanding with time. We 
interpret these structures as light echoes. 

Our formula for the density of the environment depends on 
the metallicity, here measured as the oxygen (O) abundance. 
Belley & Roy (1992) investigated the oxygen (O) abundances 
in 160 H II regions in NGC 6946 and estimated the O abun- 
dance gradient of Alog (w(O)/n(H))/Afl=-0.089 dex kpc" 1 , 
where n(O) and n(H) are the number densities of oxygen and 
hydrogen and the extrapolated log n(0)/n(H)+12=9.36 and R 
is the distance from the core of the galaxy. SNe 2002hh and 
2004et (see below) are located in —5.7 and —8 kpc from the 
center, respectively. When we adopt the solar O abundance, 
the O abundances of SNe 2002hh and 2004et are expected to 
be 8.85 and 8.65, respectively. 

Since the difference V magnitude between the peak and day 
1000 or later is -3.7, n(H) -440 cm" 3 is required if the light 
echo is associated with the SN. This hydrogen density is close 
to relatively dense gas density in such H II regions rather than 
that in interstellar space. Welch et al. (2007) suggests that 
most of the Spitzer mid-IR flux may have come from dust in 
the star formation region associated with SN 2002hh precur- 
sor. We should note that we measured the magnitude of the 
SN and the light echo. In the case of SN 2004et, which has 
the light echo and is a member of NGC 6946, the apparent ra- 
dius of the SN is —0. 1 ". If we adopted 0. 1 " radius in aperture 
photometry for SN 2002hh and we define that the inner radius 
of light echo is >0.1", -80% of the total flux is from the light 
echo as of days 1717 (our HST observations). 

Barlow et al. (2005) estimated 0.1-0.15M© dust around SN 
2002hh after 600 days with radiative transfer modelings based 
on ^/fzer/IRAC/MIPS. However, the minimum emitting ra- 
dius — 10 17 cm is too large for the emitting dust to be forming 
dust in the SN ejecta. Pozzo et al. (2006) suggest that dust 
within — 10 17 cm would be evaporated. Barlow et al. (2005) 
therefore concluded that the emitting dust must have been pre- 
existing and its origin may be from an enhanced mass-loss by 
the >10M Q progenitor. 

4.3. SN2003gd 

SN 2003gd is a Type II-P SN in NGC 628 (<i=7.2Mpc; Van 
Dyk et al. 2003) and was discovered on JD.2452803.2. Van 
Dyk et al. (2003) and Smartt et al. (2004) identified the pro- 
genitor with an F606W HST image; the progenitor is a 6-12 
Mq red supergiant. Maund & Smartt (2009) estimated the V 
and I band magnitudes of the progenitor using pre-explosion 
images from the HST and Gemini telescope archives. They 
estimated V=25.8 and 7=23.3 and an initial mass of 6-12 Mq. 

The V/?-band fading rates from day 150-300 are compara- 
ble with the radioactive decay rate as earlier mentioned. But, 
in days 300-800, the fading rates are much larger than those 
supported by the 56 Co decay power, indicating that dust for- 
mation had started from day —300 or later. Sugerman et al. 
(2006) estimated the dust mass of 1.7 xlO" 3 M Q at day 499 
(clumped dust distribution model). 

When we focus on day —500 or later, the magnitude is 
almost constant. This is due to a light echo. Sugerman 
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FIG. 5.— WHIRC view of the SNe 2002hh and 2008S field in the H- 
band. The positions of SNe 2002hh and 2008S are indicated by the white 
tick marks. A close-up image of SN 2002hh is presented in the inner box 
(2" x2"). The image was taken in a 2008 observing run (PI: M. Meixner). 

(2005; 0.225" inner and 0.375" outer radii at day 632) and 
Van Dyk et al. (2006; 0.31" at the same epoch) found the 
light echo in the F625W band. We detected an arc shaped 
light echo with -0.5" radius (5.39 xlO 19 cm at 7.2 Mpc) 
in the same band as shown in FigJT] We should note that 
the shape of the light echo is changing with time. Belley 
& Roy (1992) measured the O abundances in 130 H II re- 
gions in NGC 628 and estimated the O abundance gradient 
of Alog (w(O)/n(H))/Afl=-0.081 dex kpc" 1 , where the ex- 
trapolated log n(0)/n(H)+12=9.2 in the galaxy core. When 
we adopt R of -4.5 kpc and R{R)=232, M echo (R)-M SN (R) ^ 
10.6-2.51og n(H). Since the difference magnitude between the 
peak and day 1998 is —10.97, n(H) is —1.4 cm -3 , which cor- 
responds to a typical interstellar medium density. 

4.4. SN2004et 

SN 2004et is a Type II-P SN in NGC 6946 discovered on 
JD.2453217.5. We successfully detected the SN at all the ob- 
served bands. 

The light curves and fading rates at all bands are very simi- 
lar to those of SNe 1987 A and 1999em except for day —1000 
or later. Up to 1000 days after the initial explosion, the bright- 
ness at both optical and infrared wavelength bands was mono- 
tonically fading, but from —1200 days it starts to increase 
again, in particular in the /f-band. Our HST observations 
were taken just before and after the beginning of this bright- 
ening. Kotak et al. (2009) argued that this phenomenon is 
caused by an interaction between the SN ejecta and the sur- 
rounding ISM. On the other hand, Sugerman et al. (in prep.) 
propose a light-echo as an explanation for the re-brightening. 
Since the difference of M{V) between the plateau and day 
-1000 is -10.53, n(H) is -1.2 cm -3 . The magnitude after 
day —1000 is consistent with the light echo interpretation. 

With the bolometric light-curve Sahu et al. (2006) esti- 
mated a 56 Ni mass of —0.06 Mq. The progenitor mass is 
15-24.5 M e (Li et al. 2005; Utrobin & Chugai 2009). Fab- 
bri et al. (2011a) estimated an SN produced dust mass of 
4X10" 4 to 4xl0 -3 M Q using SED fits of combined Spitzer 
and optical/near-IR data. There are active star forming re- 
gions around SN 2004et. Like SNe 2002hh and 2008S (but a 



SN impostor, see below), the progenitor of SN 2004et would 
have formed in a metal-rich environment. 

4.5. SN2005cs 

SN 2005cs is a Type II-P SN in NGC 5194 {d=l A Mpc; 
Takats & Vinko 2006) and was discovered on JD.2453550.4. 
We note that the position measured from the F606W image 
differs significantly from the position we used for observing, 
i.e. the SIMBAD position. As a consequence, SN 2005cs was 
unfortunately on the edge of the field for the observations we 
report and was entirely missed by our NIC2 measurement in 
2008. SN 2005cs is -0.1" northeast from the bright cluster 
(Li et al. 2006). Due to its faintness, we can not detect the SN 
in the NIC2/F160W and 205W bands. 

The VRI fading rate at day 150-300 is much smaller than 
the 56 Co decay power. Pastorello et al. (2009) argued that 
such flattening could be formed by a radiation flow gener- 
ated in the warmer inner ejecta, which propagates through- 
out the transparent cooler external layers, and contributes to 
the light curve as an additional source. The lightcurves and 
magnitude decline rates beyond day —300 are comparable to 
SNe 2004et and 1987 A, suggesting that dust formation may 
have occurred. Fabbri (2011b) gave an upper limit to the 
dust mass = 3.1xlO~ 3 Mq (O-rich case) from the detection 
limit of a Gemini/Michelle 10 /im image obtained in 2007 
(PI: M.J.Barlow). 

Li et al. (2006) detected the progenitor from the HST ACS 
images, and they estimated V>25.5 and 7>24.14. Maund 
et al. (2005) and Li et al. (2006) reported the detection of 
the progenitor using the archival HST ACS images. They ar- 
gued that the progenitor is a —7-12 Mq RSG Pastorello et al. 
(2009) obtained similar results (8-15 Mq). Pastorello et al. 
(2009) also estimated that the explosion energy is — 3xl0 50 
erg (cf. — 10 51 erg in typical Type II SNe) and a 56 Ni mass of 
2.8-4.2 x 10~ 3 Mq. They measured a low expansion velocity 
(1000-2000 km s" 1 ) in optical and near-IR spectra. The flat- 
tening of the light curve in the last two epoch measurements 
suggests a light echo. Since the difference between the peak 
and day 1006 is 1 1.2 in the V-band, the ISM density is —0.4 
cm" 3 assuming that the SN has a light echo and the SN light 
after day 1006 is partially from the light echo. 

4.6. SN2006bc 

SN 2006bc is a Type II-L (Gallagher et al. 2010) or Type 
II-P (Brown et al. 2009) SN in NGC 2397 (d=\6 Mpc; Mihos 
& Bothun 1997) and was discovered on JD.2453819.4. The 
SWIFT telescope detected the SN light at an early epoch and 
the V-magnitude was —17.5 at day 15.8 (Brown et al. 2009). 
We detected the SN light in all bands. Smartt et al. (2009) 
give an upper limit of the progenitor mass of 12 M for SN 
2006bc 

SN 2006bc is surrounded by an Hll region; the optical 
spectra taken by Gemini-N/GMOS show the complex of the 
SN and H II region components (Gallagher et al. 201 1). Gal- 
lagher et al. (2011) extracted the SN component by multi- 
Gaussian fitting, and they confirmed that the Ha line profile 
is blueward shifting as the SN becomes older suggesting dust 
production. An IR-excess has been measured in the Spitzer 
IRAC bands and SED fitting by Gallagher et al. (201 1) indi- 
cates a dust mass of 9 x 10~ 4 Mq. The fading rates are much 
smaller than the typical 56 Co decay power and our estimated 
56 Ni mass of 2.7 x 10~ 2 Mq is lower than average in our sam- 
ple. A flattening of the light curve between days 400 and 500 
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FIG. 6. — (upper six panels) The temporal evolution of SN 2002hh at F606W band, (lower three panels): The temporal evolution of SN 2002hh at F814W 
band. The F606W difference images are created after PSF matching each band image. 
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FIG. 7. — The temporal evolution of SN 2()03gd light echoes at the F625W band. The echo is seen in 2004 ACS (northwestern arc structure, left panel) and in 
2007 WFPC2 observations (northern arc structure, middle panel). The radius of the echo is ~0.3" in 2004 and 0.5" in 2007. The position of the SN is indicated 
by the white stars. The difference images in right panel are created after PSF matching each band image. The radii of broken black and white are 0.3" and 0.5", 
respectively. 



suggests a short duration light echo. The estimated density of 
the matter in this light echo is ~36 cm -3 from the difference 
V-magnitude = 6.4 between day 15.8 and 694. 

4.7. SN Impostor: SN 2008S 

SN 2008S is regarded as an SN impostor in NGC 6946 and 
was discovered around JD.2454498 (Arbour & Boles 2008). 
E(B-V) by the Galaxy is estimated to be 0.34 from redden- 
ing maps of Schlegel et al. (1998). Botticella et al. (2009) 
estimated the extinction by the host galaxy Ay^l (1.13 by 
the Galaxy) by comparing their light echo model with the ob- 
served SED at 17 days. In this paper, we corrected only the 
extinction by the Galaxy with E(B-V)=0.34. Our JHK light 



curves show that SN 2008S is quite different from Type II 
SNe: for example, no plateau phase and no clear level-off 
around day ~100. The nature of the progenitor star is under 
debating. Botticella et al. (2009) argue that SN 2008S may 
have evolved from a 6-8 M Q asymptotic giant branch (AGB) 
star and formed not by Fe-core collapse but O-Ne-Mg core 
collapse by electron capture during the AGB phase. While, 
Smith et al. (2010) argue for a high mass progenitor (>20 
M Q ) for this class of objects such as SN2008S. The 56 Ni mass 
of SN 2008S has the smallest Ni mass estimated for Type II 
SNe, the value is 1.4 x 10" 3 M© (Botticella et al. 2009). 

Wesson et al. (2010) estimated the dust mass-loss rate of the 
progenitor of SN 2008S to be 5 x 10" 7 M yr" 1 and a total dust 
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TABLE 8 

Derived quantities for our targets in this study and our SEEDS sample. 



SN 


Type 


Distance 


E(B-V) 




M( 56 Ni) 


M dus , 


ISM dens. 


Light 


References for 


Name 




(Mpc) 


(MW+Host) 


(Mq) 


(xl0~ 2 M Q ) 


(Mq) 


(cm 4 ) 


echo? 


MproglM^myMdm 


1999bw 


Il-n or 
impostor 


13.7 


0.01 


? 


? 


-to- 4 


~30 


No 


(1) 


2002hh 


II-P 


5.9 


1.97 


>10 


<14 


0.1-0.15 


~400 


Yes 


(1),(2) 


2003gd 


II-P 


7.2 


0.18 


6-12 


0.8 


4xl0~ 5 -1.7xl0~ 3 


~1.4 


Yes 


(1),(3),(4) 


2004et 


II-P 


5.9 


0.41 


15-24.5 


5.9 


1.5x10^- 1.5 xl0~ 3 


~1.2 


Yes 


(1),(5),(6),(7) 


2005cs 


II-P 


7.1 


0.14 


7-12 


0.5 


<3.1xl0~ 3 


~0.4 


No 


(1),(8),(9),(10) 


2006bc 


II-P or II-L 


16 


0.52 


<12 


2.7 


9x10^ 


~36 


Yes 


(1),(11),(12) 


2007it 


II-P 


11.7 


0.13 


16-27 


9 


-lo- 4 


~19 


Yes 


(13) 


2007od 


II-P 


24.5 


0.13 


? 


0.26 


4x10^ 


~3.6 


Yes 


(14) 


1987 A 


II-P 


50 


0.19 


16-22 


7.5 


(3-5)xl0^-1.3xl0~ 3 


~1.5 


Yes 


(15) 


1999em 


II-P 


7.83 


0.10 


12-14 


1.9 






No 


(16),(17) 


2008S 


impostor 


5.9 


0.34" 


6-8 


0.14 


<0.02 


? 


Yes 


(18),(19) 



Note. — The extinction by the host galaxy is not included. The estimated dust mass is based on mid-IR data. (In SNe 2002hh and 2008S) 
Most of dust is pre-existing material. 

REFERENCES. — (1) This work;(2) Barlow et al. (2005);(3) Sugerman et al. (2006);(4) Maund & Smartt (2009);(5) Fabbri et al. (201 la) ;(6) 
Li et al. (2005);(7) Utrobin & Chugai (2009);(8) Li et al. (2006); (9) Maund et al. (2005);(10) Fabbri (2011b);(ll) Gallagher et al. (2011);(12) 
Smartt et al. (2009);(13) Andrews et al. (201 1) ;(14) Andrews et al. (2010);(15) Ercolano et al. (2007);(16) Elmhamdi et al. (2003);(17) Kozasa 
et al. (2009);(18) Wesson et al. (2010);(19) Botticella et al. (2009) 



injection into the ISM to be <0.01 M Q . The large infrared 
excess on day 17 could only be an echo, as there would be no 
way for dust in the necessary quantity to have formed in such 
a short time. The SED on day 180 could be well matched by 
models using the same shell geometry as day 17, without any 
need for new dust formation. 

5. DISCUSSION 

In Table [HJ we summarize the derived physical quantities 
of our sample. 56 Ni mass (sixth column of this table) and 
ISM density (eighth column) are from this work. The other 
quantities are from the literature listed in the last column. 

The behaviors of the SN lightcurves are affected mainly by 
(1) the explosion energy, (2) the ISM density («hX an d (3) 
the formed dust mass (Mj„ s( ) in the SNe, because the explo- 
sion energy determines the peak luminosity of the SNe, the 
surrounding ISM density/distribution and the dust affects the 
lightcurves at >1000 days and at >300 days, respectively. 
Relations between progenitor mass (M prog ) and 56 Ni mass and 
explosion energy are known for massive (>20 M Q ) Type II 
SNe (e.g., Maeda et al. 2010), however, relations among the 
56 Ni mass, M^ust, M prog , and ISM density are unknown for 
less massive Type II SNe. Through this work, we have es- 
timated the ISM density and the 56 Ni mass for our sample 
using our HST data. Using the theoretical radiative transfer 
code MOCASSIN (Ercolano et al. 2005), we have estimated 
dust masses for our sample except for SNe 1999bw, 1999em, 
and 2005cs. The dust masses listed in this table are based 
on mid-IR data and the value is as of <800 days. Here, we 



M„ 



ISM 



will examine relations among 16 Ni mass, Md ust , 
density, and E(B-V). 

5.1. Ni mass vs. M prog 

In Figj8fa), we present the relation between the 56 Ni and 
the progenitor masses M prog , . In principal, there is a relation- 
ship among the amount of 56 Ni yield in stars, the progenitors 
masses, the initial abundances, and the metallicity. The metal- 
licity of our sample is similar, we assume that the effects of 



metallicity in nucleosynthesis would be therefore negligible. 
Maeda et al. (2010) and his colleagues find relationships be- 
tween the predicted 56 Ni mass (of order 10" 2 M Q ), the SN 
progenitor masses, and the kinematic energy of the explosion 
( 1 ~ 1 00 x 1 5 1 erg s" 1 ), however they focused on Type II SNe 
evolved from >20-50 M Q progenitors (hypernova). We ex- 
pected a similar relationship between the 56 Ni and progeni- 
tor masses in Type II SNe with progenitor masses of ~ 10-25 
Mq. FigJSJa) strongly supports the previous works. As early 
mentioned, the estimated 56 Ni mass might have large uncer- 
tainties. When we exclude the data of SN 2002hh (due to 
the large uncertainty of the 56 Ni mass), M( 56 Ni) represents 
0.31 xlO" 2 x M prog .(M Q ). The 56 Ni mass derived from the 
bolometric luminosity could be an effective tool to estimate 
the SN progenitor mass. 

5.2. « H VS. E(B-V) 

FigJSfb) shows the relation between ISM density «h and 
the extinction E(B-V). We added information on SNe 2007it 
(Andrews et al. 2010) and 2007od (Andrews et al. 2011), 
which are targets in our SEEDS project. Andrews et al. esti- 
mated the nickel mass and dust mass and found the light echo 
by HST observations. The ISM densities around SNe 2007it 
and 2007od are derived by the same way applied to the other 
SNe, using the difference V-band magnitude between plateau 
and day 922 in 2007od (7.12) and day 81 1 in 2007it (8.9). For 
both objects, we assume solar metallicity. 

There is no correlation between them. This can be ex- 
plained by the different light echo geometry (shape, size, and 
inclination angle, etc.). For example, modeling of the pos- 
sible light echo geometries in SN 2004et currently suggests 
an hourglass-shaped nebula similar in size to that discovered 
around SN 1987 A, but oriented such that we are looking down 
the symmetry axis (Sugerman et al. in prep.). 

5.3. 56 M mass vs. Mdust 

FigEc) shows the relation between the 56 Ni mass and the 
dust mass Md ust . The dust mass in this diagram is the maxi- 
mum value in an early phase within l-3yr. If the correlation 
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between the 56 Ni and progenitor masses is true and we ex- 
clude SN 2002hh, where we observed the pre-existing dust 
rather than SN origin dust, there is no relation or a very weak 
relation between the dust and progenitor masses within —30 
Mq. This diagram also implies that Type II origin dust mass 
ranges from — 10" 4 upto 10" 2 Mq per SN at day <800 despite 
the progenitor mass. 



The origin of dust in high-z galaxies is a hotly-debated is- 
sue stemming from both observational and theoretical studies 
since Bertoldi et al. (2003) discovered large amounts of dust 
(-4xl0 8 Mq) in the QSO Jl 148+5251 (z=6.4, -900 Myr). 
At the present epoch in our Galaxy, Type II SNe and AGB 
stars are the main dust producers (Gehrz et al. 1989). In 
young galaxies, AGB stars are not likely to contribute sig- 
nificantly to dust production. Such low-mass stars proceed 
too slowly toward the AGB phase to produce dust within 1 
Gyr. For example, —1 Mq single stars with solar metallicity 
take — 10 Gyr to evolve into the thermally pulsing AGB (TP 
AGB). If the initial mass is —5 Mq, these stars can evolve into 
TP AGB stars and might be able to be the main dust produc- 
ers in young galaxies within —1 Gyr (Valiante et al. 2009). 
However, massive stars can evolve into Type II SNe in times- 
pans under 20 Myr. Theoretical ISM dust models by Dwek 
& Cherchneff (2010) predict that an average 20 Mq (initial 
mass) SN has to make at least —0.1-1 Mq of dust in order to 
be a viable source for the dust found in the QSO J 1 148+5251. 
However, the 56 Ni mass-Mj„ s/ diagram does not follow the 
prediction by Dwek & Cherchneff (2010). 



6. SUMMARY 

We performed BVRIJHK band photometry of 6 Type II 
SNe at late epochs (days >500) using HST and ground 
based telescopes to investigate their natures and relations be- 
tween dust mass, 56 Ni, ISM density, and progenitor mass. 
Most of our HST observations are successful in detecting 
the light from the SNe alone and in measuring magnitudes 
with less contamination from nearby stars. Combining our 
data with previously published data, we showed lightcurves 
at VRIJHK-bands and estimate the decline magnitude rates 
at each band at 4 different phases. These lightcurves and 
other data indicate that dust is forming in our targets from 
day —300-400, supporting SN dust formation theory. We es- 
timated the ISM or circumstellar density around SNe 2002hh, 
2003gd, 2004et, and 2006bc. The light echo density around 
SN 2002hh is higher than a typical ISM density. The ISM 
density around the other SNe is a 1-70 cm" 3 . We estimated 
56 Ni masses (0.5-14 x 10~ 2 Mq) by comparing their bolomet- 
ric luminosity with SN 1987 A and we find that it correlates 
with progenitor mass. This relation supports the work by 
Maeda et al. (2010), who focus on Type II hyper SNe (>20 
Mq). The dust mass does not appear to be correlated with 
56 Ni mass among the 7 SNe. 
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